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1 Introduction 



The currently available experimental data on deep inelastic structure functions covers a large kine- 
matical regime with high precision measurements. This provides an interesting challenge for theo- 
retical physics when it comes to describing this data in the low-Q^ domain. At sufficiently high ^ 
a perturbative approach is justified and within the classical DGLAP approach next-to-leading-order 
(NLO) fits are still the state-of-the-art. Although some authors have been using NLO DGLAP 
down to rather low values of , it becomes more and more obvious that such an approach is not 
sufficient for a precise description of deep inelastic structure function data. Missing are not primar- 
ily the sub-leading terms in powers of Og, but corrections which are proportional to the reciprocal 
value of the photon virtuality Q^, viz. higher-twist terms. This statement gets some support from 
the experimental observation that there is a rather large number of events in which a longitudinally 
polarized photon scatters off a proton diffractively and results in a vector meson. Since this process 
has been shown to be calculable perturbatively and to belong to higher-twist one should 

expect a sizeable amount of higher-twist contributions to deep inelastic structure functions in the 
same kinematical regime. However, it has been shown recently that the situation is more compli- 
cated Q: one gets contributions with different signs which might lead to strong cancellations and 
an estimate of the amount of inclusive higher-twist contributions based exclusively on diffractive 
vector meson production is not sufficiently reliable. Therefore a systematic study of higher-twist 
is not only worthwhile but necessary. This is even more true, since at least some of the twist-four 
contributions, after removing the explicit 1/Q^-dependence, are known to evolve faster than the 
leading-twist contributions, and to increase faster with decreasing xb- 

Most computations based on twist-four operators for the inclusive process have been published 
in the early 1980s. S. Gottlieb classified all twist-four operators and computed some of their 
anomalous dimensions. The anomalous dimensions of four-quark operators and of flavor octet, 
three-body operators were obtained by M. Okawa |10|, and the Wilson coefficient functions of four- 



quark operators by S.P. Luttrell et.al. Calculations based on the complete set of twist-four 

operators introduced in Q are difficult since this set is over-determined: different operators are 
related by the equations of motion. This deficiency was remedied by R.L. Jaffe and M. Soldate 
|l2| , p!^ who introduced a special basis, which they called canonical. They showed that this basis 
is adequate for the necessary twist-four perturbative calculations and calculated the operators' 
tree-level coefficient functions. 

These approaches are all based on the operator product expansion (OPE) and make use of the fact 
that the hadronic tensor W"'^ factorizes into two parts, viz. coefficient functions Cj and parton 
correlation functions fi which depend on a set of light-cone momenta {x}: 

T-2 

^"''= E {$) ' E (1-1) 

T even 

where r corresponds to the twist of the quantities and we have dropped details like the dependence 
on logarithms of Q^. On the other hand, R.K.Ellis, W. Furmanski and R. Petronzio (EFP) Q 
followed a more intuitive approach and computed 1/Q'^-corrections to the conventional leading- 
twist Born diagram starting from leading-order Feynman diagrams. Using special techniques to 
classify the spinor, color and momentum structure, the authors could perform the factorization of 
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quark-coefficient functions and operators and reduced the number of contributing terms witli the 
help of equations of motion. 

The computations in ^ can coarsely be divided into two steps. In the first step the authors 
considered the relevant diagrams with four f-channel partons at the Born level and performed the 
factorization of the coefficient functions and matrix elements of the relevant operators. In detail, 
these operators are (we are ignoring terms which involve a matrix 75 )0: 

^^(A) = ^{p\T {mYH>^n)} \p) , (1.2) 
B'^ivA) = \{p\T{mY^^ir]nm^n)} \p) , (1.3) 
Cf^^'ito, r], A) = ^{p\T {i^{0)YB^{i^n)D''{'nn)iP{Xn)} \p) (1.4) 



and 



Qico, r/, A) = ^/(P|T { [mfit^'H^n)] [^Pir^n)f^t''^PiXn)] } \p) . (1.5) 

In the second step, EFP used equations of motion in order to express the contributions of the two- 
quark operators ([1.2|)-(p~4|) in terms of linearly independent operators. For this they studied both 
the longitudinal basis (which corresponds to the canonical basis in |jl2|) and introduced a new one 
which they called the transverse basis. While within the longitudinal basis the expressions turned 
out to be rather complicated, in the transverse basis the twist-four contributions of the two-quark 
operator to the transverse and longitudinal structure functions Ft and Fl can be expressed in 
terms of 

A2Ti(x) = ^e'^-(p|T {v;(0)7^^7.d^^'D^K0)d'^-'D,,(An)V^(An)} \p) (1.6) 



and 



A2T2(x2,xi) = ^^e^''(^^-^^)e^^^Mp|T{v5(0)7;.^7.rf""'D.'(r/ri)d'^'^'D^,(^™)V'(An)} \p) , 

(1.7) 



such that the twist-four contributions to -Ft and Fl are (we are neglecting the proton mass) 



A' 



2 

TpT=4 _ 

A' 



4Ti(a;B)-XB /dxidx2 T2(x2,xiJ 

X2 - Xi 



(1.8) 
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= A-^TiM (1.9) 



plus terms involving the four-quark operator (1.5). 



EFP's results have been used as a starting point for many other studies. For example, introducing 



a quantity named special propagator, J. Qiu [14| could reproduce the results obtained in |j^ in a very 

^Cf. Eqs. (2.45)-(2.50) and (5.23)-(5.24) in The calculations are performed in the axial gauge with the gauge 
vector n^. Since we are interested in the general structure of the results, we postpone the definition of various 
quantities until later sections. 
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elegant way. A. P. Bukhvostov et.al. ||T^ define a special class of operators which they call quasi- 
partonic, formulate the evolution equations for these operators and compute the integral kernels 
for all possible 2^2 transitions in the one-loop approximation. 

All these computations of twist-four quantities are based on quark-operators. In the small-XB 
domain, on the other hand, we know that in the leading-twist case the gluon-operator is dominating 
since it leads to the maximum power in I/xb- For higher-twist, one expects that in the small-XB 
domain it is, again, the gluon-operators which lead to the main contributions. A systematic analysis 
of twist-four gluon-operators is missing. Such an analysis could, again, be based on the OPE. For 
twist-four one has to take into account contributions of two-gluon operators, 

trF^i"D^^..D'^"-iF%g^,^^ + perm., (1.10) 

as well as of four-gluon operators. 





7?Mi 
• a • ■ 




. . F^' 


trF'^i". 




. . r a ■ 


. . F^ 






pH 


. F^' 



Here the dots denote products of covariant derivatives. In a future analysis one will have to compute 
the corresponding coefficient functions and the anomalous dimensions. This is a demanding task 
since all these twist-four gluon-operators will mix under renormalization. 

Alternatively, one can make use of the more intuitive approach of EFP in order to find expressions 
for the contributions of twist-four gluon-operators to inclusive DIS structure functions. This is 
what we want to discuss in this article. We base our computations on the diagrams in Fig. || 
and consider amplitudes with two, three and four t-channel gluons. The method we are going to 
describe will lead to the factorization of gluonic coefficient functions and matrix elements of gluon- 
operators between proton states. For our analysis we are neglecting the target-mass corrections 
and we do not start from the classical way of defining twist as the difference of canonical dimension 
and spin; instead we will search for l/Q^ corrections to the familiar leading-twist result. We will 
use a dimensional analysis to show that diagrams with more than four t-channel gluons lead to 
contributions of, at least, twist-six. We are working in the axial gauge, and we can show that, 
within this gauge, it is sufficient to restrict oneself to the diagrams in Fig.|l|. In this paper we do 
not compute explicit expressions for the upper quark-loop parts in the diagrams of Fig. |l|. Instead 
we derive some general expressions for the corresponding amplitudes and make use only of the 
color-, Lorentz- and momentum structure of these quark- loops. 

In comparison to quarkonic operators considered by EFP, the case of gluonic operators is slightly 
more complicated, and we will carry out only the very first steps of the corresponding analysis. 
Since the photon couples to the four-gluon state only through the quark loop, already the "Born" 
approximation will contain contributions proportional to logQ^ (i-e., UV-divergences) . In this first 
part of the analysis, we will ignore logarithms and we will not address the computation of the quark 
loops. The emphasis of the present paper lies on the analysis of the color and Lorentz structure and 
on the factorization of the quark loop and the proton matrix elements. We will not yet consider the 
moments of the structure functions and also the definition of an operator basis will be postponed 
to a future paper. 

For our discussion it is useful to clarify the notion of leading-order (LO) and next-to-leading order. 
The diagrams in Figs.|^ and c are suppressed by factors g and g'^, resp., compared to the diagram 
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in Fig. Within our derivations we will factor these constants out and absorb them within the 
expressions for the lower blobs. This way we can treat all three quark-loop expressions as if they 
were of the same order in Os- On the other hand, computing these quark-loop expressions one will 
get terms proportional to logarithms in and terms which are just constants. Terms which are 
proportional to logarithms correspond to LO quark-coefficient functions convoluted by Altarelli- 
Parisi splitting functions, while we denote the constants as NLO gluon-coefficient functions. For 
this first analysis we restrict ourselves to quark-loop diagrams which are one-particle irreducible 
(IPX). The three- and four-gluon vertices in Fig. ^ have, therefore, to be considered as if they were 
absorbed by the lower blobs and are part of the evolution. 

For our computations it is convenient to use the axial gauge, i.e. we are introducing a light-like 
vector n'^ such that 

n^A^^ = 0. (1.12) 

We are using the same vector n within the decomposition of momentum vectors into Sudakov 
variables: 

/cf = Xip'^ + ai'n'^ + A;^rp , = = p ■ ki^q^ = n ■ ki^i = , p ■ n = 1 . (1-13) 

Here p represents the momentum of the target hadron, which we assume to be massless throughout 
this paper. For the color factors we have chosen the following normalization of the Gell-Mann 
matrices i": 

tnH^ = ]^5''\ (1.14) 

and we set the number of colors N^. = 3. 

The outline of this paper is as follows. In section |2| we define our basic quantities. In particular 
we distinguish between the upper quark-loop parts in Fig.||, R, and the lower blobs, F. The 
expressions for the quark-loops will then be expanded about the longitudinal components of the 
gluon-momenta they depend on. For the derivatives which occur in these expansions, we compute 
some Ward identities in section |^. In section ^ we show how, in the axial gauge, the ^a-operators 
which occur in the matrix elements for the lower blobs in Fig.|l| can be expressed in terms of Fa^- 
operators. Section || is devoted to the dimensional analysis of the quantities under consideration. 
This analysis leads to the justification that we can restrict ourselves to the diagrams in Fig. || 
for a consistent analysis of gluonic twist-four contributions. Since the diagrams in Fig.|l| contain 
also leading-twist contributions and contributions of twist larger than four, we will provide, also 
in section ^, a method for projecting out the twist-four parts. Using the results of sections § 
through ^, the twist-four contributions we are interested in can be written in a completely gauge 
invariant form, which is factorized into coefficient functions and matrix elements between proton 
states. We give these results in section |6|, viz. Eqs. (|6.5|), ( |6.1[)| ) and ( [6.251 ). Subsection 3^ contains 



a comparison between our results and those obtained by EFP. Finally, in the conclusions, we 
discuss the steps which have to be completed before our results can be used in a phenomenological 
analysis. 



2 Factorization 



As outlined in the introduction it is our aim to specify twist-four gluon operators, write them 
in covariant form and perform the factorization of coefficient functions and matrix elements of 
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(a) (b) (c) 

FIG. 1: For a systematic analysis of the twist-four contributions of gluonic operators it is sufficient to consider 
the contributions of the three diagrams shown here. The quark-loops symbolize the coupling of the 
gluons in all possible ways and the propagators of the gluon-lines belong to the lower blobs. We 
take all gluon momenta flowing upwards, i.e. into the quark-loops. 



these operators between proton states. We will argue that it is sufficient to consider the diagrams 
in Fig. 1^. Here the upper parts represent the quark-loops with two, three and four t-channel 
gluons. They are amputated in the gluon lines, and we denote them by R'j^^{ki), R'j^J'^{ki, k2) and 
R'^^^Xy{ki, k2, ks), resp. These quark- loops will contribute to both NLO twist-four gluon-coefficient 
functions and terms which are convolutions of LO quark-coefficient functions with Altarelli-Parisi 
splitting functions. It is our aim to derive general results for the twist-four contributions of the 
diagrams in Fig. || without using explicit expressions for these quark- loops. All we need are some 
general properties of the functions R concerning the Lorentz-, color- and momentum-structure. 

The lower parts in Fig. |l] are denoted by r^^(A;i), ^abc{ki, k2) and ^^cd ^s)) resp., and with 

these conventions the contributions of the diagrams in Fig. |l] are given by 

j^,Rtu{ki)TZ{kr) , (2.1) 

^2(xB,Q') = I^^R;tih,k2Kb:{h,k2), (2.2) 

Hs{x^,Q') = / (2?)4 (1^)4 (t?)4 <t k2, k.K-^^Jik,, k2, fcs) , (2.3) 

where we have dropped the photon indices. As we will show in section ^ expressions ( ^.lD -( ^3|) 
contain leading-twist and twist-four contributions of gluonic operators as well as some parts with 
twist larger than four. The twist-four parts can be projected out if one performs a dimensional 
analysis. Here, the convenience of the axial gauge becomes most explicit. 

The lower blobs in Fig.|l] can be written as matrix elements of gluon-operators between proton 
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FIG. 2: Sample diagrams which contain three- and four-gluon vertices. 



states: 



abc {kl,k2) 



d zi e 



iki-zi 



{p\T{Aim^aizi)}\p) , 



(2.4) 
(2.5) 



d zi d ^2 d 23 e 



ifel-zi i(fc2-fci)-Z2 J{k3-k2)-Z3 



{p\T{A^,m'Az3)AUz2)A^aizi)}\p) . 

(2.6) 



Here T denotes the time ordered product, and expressions (p.4|)-(^^) contain the propagators of 
the gluon hnes in Fig. ||. 

In order to perform the factorization procedure with expressions (2.1)-(2.3) it is convenient to ex- 
pand the expressions for the quark-loops about the longitudinal components of the gluon-momenta 
(cf.dTTsD): 

a2 nab , 



, , dR±ixip) 



1 d^RZixip) , 
{ki - xipY + - (ki - xipYiki - xip)^ + ■■■ , 



2 dkldkt 



R%l[k^, k2) = <^(xlp, X2P) + ■ 



dRfauixiP,X2p) 



i=l 



(ki - Xip)^ + 



Rf^Aki, k2, ks) = ^2P, X3P) + ■ 



abed 



(2.7) 
(2.8) 
(2.9) 



We will justify in section ^ that we can neglect additional terms in these series if we are interested 
only in terms up to twist-four. 

Before we give the results of the factorization procedure, we compute some Ward identities for the 
derivatives in Eqs. ( p.7D and ( |2.8| ), show how the expressions for the lower blobs F can be expressed 
with the help of gluonic i^^'^-operators, and perform a dimensional analysis of the expressions 
(|2.l|)-(l2.3|) in the next three sections. 



3 Ward Identities 



In order to find relations for the derivatives of the quark- loops in Eqs. ( ^.7] ) and (2^) we will 
make use of Ward identities. Although these identities can be derived diagrammatically at the 
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one-loop level, this is a rather tedious task. We, therefore, make a digression into the path-integral 
formulation of QCD and derive these Ward identities in the most general way. For this, we first 
derive a relation for the generating functional of perturbatively calculated vertex functions V in 
the next subsection. This relation, viz. Eq. (|3.9D , will then be used to derive the desired equations 
in the following two subsections. 

We should stress that the Ward identities which will be derived in this section are valid for IPI 
vertex functions. They can, therefore, not be applied to expressions for quark-loops which involve 
three- and four-gluon vertices like those in Fig.|2|. 



3.1 Generating Functional 

The generating functional for Green's functions in the axial gauge can be written as |1(:, 17 

W[J^,K'']=J^ J^im" exp!^i y"d^x [^cD - A'^n-^° + A;jJ^ + A'^i^'^jj . (3.1) 

Here 3sf is some normalization factor, and K"^ are sources, the term X'^n-A'^ has been introduced 
in order to define the axial gauge. The generating functional must be invariant under gauge 
transformations, which for infinitesimal uj{x) can be written as 

A;^A; + (A)a;^ = A; + d^u;- - gr'^^A^J . (3.2) 

Inserting (|3.2|) into (|3.1|) , substituting fields that occur as factors within the integral by functional 
derivatives and making use of the arbitrariness of w, one gets 

= J^f{J) (n'^T^^ - Ja^)) WIJ^^K'^] , (3.3) 

where the covariant derivative is given by 

Df{J) = 5-'d,-gr'^^. (3.4) 

For our analysis we need the generating functional V for vertex functions. This can be obtained 
with the help of the following procedure. First, define the functional V via Legendre transformation 
of the generating functional of the connected Green's functions, Z, 

V[AlA''] = Z[r^,K'']- jd'x [A'l.Ji: + KaA''] , (3.5) 

where Z is defined as 

Z[r^,K''] = -ilnW[r^,K'']. (3.6) 
With these definitions we get from (|3.3| ) 

D^^\A)L^^A, + ^] =0. (3.7) 
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Now, V can be obtained via subtraction [|l^,p^: 



and the Ward identity (3.7) finally becomes 



, 6V 



0. 



(3.5 



(3.9) 



3.2 Relations Between Two- and Three-Gluon Quark-Loops 



We are now deriving an equation which relates the derivative of the two-gluon quark-loop with 



the three-gluon quark- loop R^^i^- The relation between the perturbatively computed quark- loops 
and the generating functional V is as follows (see Fig. || for our conventions on the gluon momenta: 
the nth momentum 4- vector kn, i.e. the nth argument of the i?- functions, is the sum of momenta 
of the first n gluons): 



J d^x d^xi d^X2 e*'=i-^ie'('=2-'^i)-^e''=-^2 



6A^{xi)6A-{x2) 
6A'iixi)6Alix)6A-ix2) 



{2TT)^6{h + k)R'±{h) , (3.10) 



(3.11) 



/d^xd^a;id^a;2 d^^s e''=i-^^e^(^2-fci) Xgi(fc3-fe2)-x2gifc.x3 

J OA'a 



(X1)M^(X)M,^(X2)M2(X3) 

--{27rf6{ks + k)R'tt{k^,k2,k^)^ 



(3.12) 

Relabelling Lorentz- and color indices in Eq. ( |3.9| ) and differentiating the equation by Aa{xi) and 
A'^{x2) one gets 



6Aii{xi)SA-ix)6A-{x2) 



g { r'"' 5{x - xi) 



6A'^,{x)6A^^{x2) 



f^' 6{X - X2) 



5A>'a{xi)5Al,{x) 



(3.13) 



Now, multiplying Eq. ( 3.13| ) by appropriate exponentials, integrating over x, xi and X2 and keeping 
Eqs. ( p. 10 ) and (3.11) in mind one ends up with 



-ik'^Rf,',{ki,k + k^)=g[r^''Rlt{k + k^)-f''''R%{k^)] 



(3.14) 



The two gluons at the lower end of the quark loop must be in the color singlet. We can, therefore, 
define 

Rf^{k)^5'''>R,,{k). (3.15) 
With this definition we can take the limit A; — > in Eq. ( p. 14]) and finally get the desired result: 

(3.16) 



■^9r''^-^^ = RtAkiM 



dkl 
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3.3 Relations Involving Four-Gluon Quark-Loops 



Relations between derivatives of the three-gluon quark-loop -R^^'L the four-gluon quark-loop 
^liaXu obtained in a completely analogous way. If we differentiate (|3.9| ) three times, multiply 

the resulting equation by suitable exponentials and perform the Xj-integrations, we get, similar to 

-ik'^R'jii.ih, h + k,k2 + k)= g{r'''R'Xiki + k,k2 + k) (3.17) 

- f'^'R^itikiM + k)- f"''' RltikiM)} . 

Before we can take the limit /c — > we have to explore the color structure of R^^^y First, consider 
the one-loop case. Since -R^JtL is the sum of all diagrams with three external gluons, we can, in 
the one-loop case, always group two diagrams together. For each diagram there is another diagram 
which gives the same momentum dependence (up to the sign) but different color structure (see 
Fig.|^ for illustration). Since the gluon-system must have even C-parity this color structure is 
always of the formQ 

ti{etH^) - tr(t^tV) = i i (3.18) 
for each pair of diagrams. We, therefore, define 

<A^(fcl, k2) ^ gr'^RMh, k2) , (3.19) 

where the factor g has been factored out for convenience. We expect that Eq. ( 3.19| ) reflects the 
general ansatz for the color structure of a three-gluon system in the case of even C-parity and is 
valid beyond the one-loop level. If we, now, insert ( P^ into (|3T7|) and make use of the Jacobi 
identity, 

jabc' jc'cd jacc' jc'bd _j_ jadc' ^c'bc g 2Q^ 

we can take the limit A: — > and obtain 

■2 \ fahc' fc'cd(^^IJ.>'v{ki,k2) ^acc' /.c'fed ^-^M^i^ (^1 ) ^2) 1 jyahcdd. . r \ /'Qr,1^ 

15 I J J + J J ■Q^ j = R^,a\v[kl,kl,k2) . (3.21) 



For our computations we need relations for the derivatives of R'^^x^ with respect to ki or k2- Since 
Eq. ( 3.21| ) contains these derivatives only as a sum, we have to project with suitable color tensors. 



Now, we realize that there is more than one color tensor which projects out the desired quantity 
in Eq. ( |3.21| ). For example both 

pabcd ^ ^ac^bd (3 22) 

and 

pabcd ^ ^acc'^c'bd (3 23) 
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In the case of odd C-parity the coupling of the gluon to quarks and anti-quarks involves a sign change. In this 



case, one therefore gets a factor tv{t'^t''t'') +tr(t''t'T) = ^d"*"" 
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FIG. 3: Two sample diagrams which enter into i?"^'^^ at leading order in as- With an adequate choice of the 
loop integration variable the traces over Dirac-matrices give the same (up to the sign) momentum 
dependence for both diagrams. The first diagram contributes with a color factor tif^t^f^ while the 
second one contributes with a color factor tr t'^t^'t'^. 



project out dRfj,xu/dkf . This fact will lead to some important restrictions for the four-gluon quark- 
loopi?^tt (Eq-O below). 

For the remainder of this section we will restrict ourselves to the one-loop approximation of the 
quark loop. In this case we can use the same arguments that led to the definition ( |3.19D . For four 
gluons we define the quantity (cf. ( |3.1^ )) 

= _lfabc'jc'cd _^ ^^abc' ^c'cd ^ \§ab^cd _ (o ^/^^ 

4 4 6 

We realize that there are three linearly independent tensors d, which can be obtained by permutation 
of color indices in Eq. ( |3.24| ) and, therefore, write the color structure of R'J^^^u 

Rf:iiku k2, k,) ^ {d^'^'^Rl^x^ik,, k2, k,) + d'^^'^^Rl^xuiki, k2, k,) + ^-'"^/^^^..(fei, k2, k,)} . 

(3.25) 

Each diagram that contributes to R^^^x^, at the one-loop level enters into one of the functions R^^^xu • 

Inserting ( |3.25| ) into ( p. 21 ) and projecting with the tensors p^^'^'^ and P!^^'^'^ we end up with two 
equations, viz. 

^^M^gi^ ^ -^«i.A.(fci,^i,fc2) + ^^R%xAkiMM) + ^^Rl^xAkiMM) , (3.26) 



20i = -^^RUuihM, h) + ^^RUAkuh, k,) + ^-^RUAkuku k,) . (3.27) 



These two equations can simultaneously be fulfilled only if the following equation hold^: 

RlaxAki, kiM) + Rl^xuiki, ki,k2) + Rl^xuiki, ki,k2) = 0. (3.28) 

Note that this equation holds only true if one of the gluon- momenta vanishes (see Fig.|| for our 
conventions on the gluon- momenta) . Inserting the latter equation into Eq. ( p. 261) we end up with 
the desired results: 



Qp<^b<^ (k k ] 1 

i ^ g^.'' = --gr'^RUAkiMM) = --gr"'RUauiki,k2,k2) (3.29) 



^There are more tensors than those in Eqs. ( |3.22| ) and (3.23) which project out dR^w /dki in (3.21). Using these 
gives a consistent picture but no new results. 
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and 



dR±l(ki,k2) 1.1 1 , 



(3.30) 



Here we have also projected out the derivative of R'f^J'i, with respect of k2 in Eq. (|3^ ) and we have 



used the fact that we have some freedom in choosing Lorentz- and color-indices and momentum 
variables in defining the four-gluon quark-loop in Eq. ( p.l2| ) . R^^f^Xf ™ust be invariant under the 
simultaneous exchange of momentum variables, Lorentz- and color indices. This leads to the 
possibility to express the derivatives of -R^^fi/ ™ ^qs- ( 3.29 ) and ( 3.30 ) either by Rj^^xu °^ RfiXau- 
We are discussing this symmetry in more detail in section 

Finally, differentiating Eq. ( |3.16| ) by k^ and using (3.29) and ( 3.3C| ) we get 



dRfAh) 



dkldkl 



b"" {RUuikuh, fci) + Rl^xuih, fci, k,)} 



(3.31) 



where we have also made use of ( |3.28| ). 



4 Gluonic Operators 

As outlined in sectio n ^, the lower blobs in Fig. ^ can be written as matrix-elements of gluonic 
operators (Eqs. ( |2.4| )-( |2!q )). Since it is our aim to find expressions for the gluonic twist-four con- 
tributions which are explicitly gauge invariant, we have to find relations in which all operators 
consist of gluonic operators Fa'^ and covariant derivatives D^^j^. In this section we show how, in the 
axial gauge, Aa-operators within matrix elements can be substituted by F^^'-operators. For this 
we make use of the relation 

n^F^ = (n•a)A^ (4.1) 

Let us, first, define^ 

jd'zd'z'^^-^e-'^'-^\p\T {Al{z')A>:^{z)] \p) ^ {2,r)H{k - k')GZ{k) , (4.2) 

and examine the expression 

y"dSdVe''=-^e-''^'-"'(p[T{(n-a')Ab(/)(n-a)A^(z)} \p) = 

= -in ■ kin ■ k' jd^z d^z' e'^-'e-'^'-'' {p\T {A'^{z')A^{z)} \p) = 
= n-kn- k'{2TTf5{k - k')G^^{k) , (4.3) 
®We are not writing the p-dependence of G^^ explicitly. 
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where we have used partial integration. If we perform the integration over k' in Eq. ( |4.3| ), keeping 
(|4.1| ) in mind and inserting x = n ■ k, we end up with 

Jd^ze^^-^iplTiA^^mi^iz)} \p) = ^ y"dSe''=-n„n^(p|T {f^ {0)F:^^ (z)} \p) , (4.4) 

which is the desired result. In complete analogy, one can find relations for matrix-elements with 
more than two gluon-operators. 



5 Dimensional Analysis 

Before we present the results of the factorization procedure we, first, perform the dimensional 
analysis in this section, in order to justify that we have restricted ourselves to diagrams with up to 
four t-channel gluons (Fig.|l]) and to those terms in the expansions of the quark- loop expressions 



R which are given in Eqs. (|2.7D-(2.9). By this, we also derive the projectors which are needed to 



obtain contributions of definite twist from the expressions under consideration. 

To understand the reasoning of the following discussion, it will be sufficient to anticipate that after 
factorization the upper and lower parts of the diagrams in Fig. || are connected by integrals over 
momentum fractions Xj, which we can write as 

Hi = y"dxii?^,(xip)2G^"(xi), (5.1) 

H2 = jdXldX2 R^au{xip, X2P) {Xl, X2) , (5.2) 

= y dXidX2dX3 R^aXuixiP, X2P, X3p) ^C"^" {xi, X2, X3) . (5.3) 

It will also be sufficient to know that the functions G are matrix elements of gluon operators which 
may contain derivatives. The number of Lorentz-indices corresponds to the number of gluon-lines 
plus the number of derivatives. For convenience we are using the above relations already in this 
section. The precise definition of the functions G will be given below in Eqs. (U), ( |6ll| ) and (|636D 



and the functions H^, which are not the same as the functions Hi in ( |2.1[ )-( p^ , will be discussed 



We expect an infrared cut-off to enter into the expressions for the quark-loops only within loga- 
rithms since infrared and mass singularities are logarithmic at most. Moreover, renormalizability 
guarantees that one-loop diagrams with two external photons and any number of gluons do not con- 
tain UV divergences. Therefore, after factorization, the only dimensional quantity the quark-loops 
R can depend on is Q^. On the other hand, the only dimensional quantity the matrix-elements 
G can depend on is A. If we, therefore, determine the dimension of a particular function R with 
Lorentz indices projected out, we can relate this dimension to the powers of this function 

must be proportional to. 

A simple analysis shows that the dimension of the quark-loops is given by the expression 

dimR = 2- B , (5.4) 
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where B is the number of gluon-hnes plus the number of derivatives: 

B = UG + riY) . (5.5) 
With the help of (|riD-(|5j) we conclude that 

dim 4G = dim 3G + 1 = dim 2G + 2 . (5.6) 

5.1 Two Gluons 

In order to split up Hi, ( |5.lD , into leading- and next-to-leading twist expressions, we, now, explore 
the tensor structure of 2^^^^ (which is given by Eq. ( |6.6| ), below)|^ 

^G^^ixi) = d^^'Aiixi) + n'^n^^2(xi) . (5.7) 

Since we are working in the axial gauge, i.e. n ■ A = 0, there can be no terms proportional p^p'^ in 
(|57|) (cf. (U)). The tensor 

has been introduced since it serves as a projector, such that 

Ai = ]^d^,2G^\ (5.9) 

A2=p^Pu2G^'' . (5.10) 



Now, since dimp = 1, dimd^'^ = and dimn = —1, we conclude with the help of ( ^.61) that 
dim ^2 = dim ^1 -|-2. It follows that in Eq. ( ^.7| ) the term proportional to Ai gives the leading-twist 
contribution, while the term proportional A2 contributes to twist-four, i.e. in (|5.1| ) the corresponding 
expression is proportional to Is? /Q"^ . 



5.2 Three Gluons 



The twist analysis of ( |5.2D ca n be performed im complete analogy. Exploring the tensor structure 
of aG'^""'^ (defined in Eq. ( |6.1lD , below) one getsQ 



Here, we have made use of the fact that 



0. 



(5.11) 



(5.12) 



As mentioned before we postpone the derivation of some of the results we are using here to the 
next section. Eq. ( ^.12 ) follows since the indices of sG^'^'^ belong either to terms like A^ or like 
where the tensor u;^, is introduced in (|6.7D. It subtracts the longitudinal component of the 



^Additionally, terms proportional to UaPise^""^ (in Eq. (5/7)) or naS^""" (in Eq. ( 5.11 )) would be present in the 
case of polarized scattering or if a virtual Z boson scatters off the proton. 
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momentum vector it is applied to and UaUJ^' = 0- The various terms in ( 5.11 ) can be projected out 
in the following way 

^N^Vp-sG^"" etc., (5.13) 
B2=P^.PuPa3G^'''^' , (5.14) 

and a dimensional analysis shows that the terms proportional Bl contribute to twist-four, while 
the term proportional B2 contributes to twist-six. 

5.3 Four Gluons 



The tensor structure of i^G^"^'^ (see Eqs. ( 6.26| )-( |6.28D below) is given by 



^Qt^aXu ^ f^t^y^oXfji ^ pg^^^ ^ d^"'n''n^Cl + perm. + n^vJ^n^'n^C^i , (5.15) 

with 

C| = d^.<xA4G'^"^" etc., (5.16) 
G2' = ^Vp-Pa4G^"^" etc., (5.17) 
C3=P^P.PaPA4G^"^^ (5.18) 

where 

d^vaX = o {"^d^iydfjx — d^adux — d^xdua) ■ (5.19) 
o 

In this case, the functions CJ give the twist-four and the functions the twist-six contributions, 
while C3 belongs to twist-eight. 

Terms with more than four t-channel gluons and/or derivatives can be analysed in the same way. 
The term with lowest dimension of the corresponding five-gluon contribution r^G^"^'"^ would be 
proportional d^'^d'^^n^ + (perm.) and, therefore, belongs to twist-six. Increasing the number of 
gluon-lines or derivatives also increases the twist. We can, therefore, justify that we have restricted 
ourselves to the diagrams in Fig.|l| and to the terms in Eqs. (2.7)-(|2.9|). 



6 Results of the Factorization Procedure 



Before we present the results of the factorization procedure, let us summarize the preceding results. 
Our computations are based on the diagrams in Fig.|l|, where the expressions for the quark- loops 
R are expanded about the longitudinal components of the gluon-momenta ( p.7D -( ^^ . For the 
derivatives in ( p. 7] ) and ( p.g| ) we have computed some Ward identities in section ^, viz. ( |3.16| ), 
( 3.29| ), ( 3.30 ) and ( 3.31| ). For the derivation of these relations we had to explore the color structure 
of the quark-loop expressions, (p.l5| ), ( p. 19 ) and ( |3.25| ), which also led to an important relation 
for the four-gluon quark-loop, (|3.28D , if one of the gluon-momenta vanishes. In section Q we have 
developed a method which is useful in order to replace, within matrix elements, ^d-operators by 



14 



F-operators, and in section ^ we performed the dimensional analysis of the expressions under 
consideration. This yielded a scheme for projecting out contributions with definite twist. 

As outlined before, it is our aim to write the gluonic twist-four contributions to the inclusive deep- 
inelastic structure functions in an explicitly gauge invariant form. In order to get these expressions 
it is, now, necessary to reorder the terms which enter into the sum of the three diagrams in Fig. || 
after expanding the quark loops about the longitudinal components of the gluon-momenta as in 
Eqs. (p.7|)-(pl9|). In particular, we write the sum of these contributions as 



Hi + H2 + Hs = Hi + H2 + H3 + 0[ 



(6.1) 



Here Hi consists only of the zeroth-order term of the Taylor expansion of the two-gluon quark-loop 



in Eq. (2.7) inserted in (2.1), while H2 is the sum of the first-order two-gluon quark-loop (from 
Eq. ( |2.7| )) and the zeroth-order three- gluon quark-loop (from Eq. ( |2.S| )) inserted in Eqs. (2.1) and 
(|2.2| ). consists of the sum of expressions ( |2.lD through (2^) where only the zeroth-order term of 
the four-gluon quark-loop (p.9|), the first-order three-gluon quark- loop (|2.8[) and the second-order 



two-gluon quark-loop (2/7) are taken into account. With these conventions the steps needed to 
write Hi and H2 in an explicitly gauge invariant form are straightforward, while the situation is 
more peculiar in the case of H3. We review the results in the next subsections. 



6.1 Two and Three Gluons 



Inserting the zeroth-order term in Eq. (2/7) into ( |2.1[) , we realize that the expression for the quark- 
loop depends only on the longitudinal component of the four-momentum ki, and we can simplify 
the integrations over ki and z if we introduce the expression (cf. ( |1.13D ) ||l| 



and make use of the relation 

d^k 



(2vr) 



^e'''' '6{xi -ki-n) 



^e'^-i5{4)(^_An). 
27r 



The result is 



dA 



Hi = IdxiR'liixip) I :^e'^'^Hp\T{AmA^ai>^n)}\p). 



(6.2) 



(6.3) 



(6.4) 



With the help of the results of section we can substitute the ^-operators by -F-operators, and 



with the help of Eqs. (|5.7| ) and ( 5.10 ) we can project out the twist-four contribution: 

Hi\^^^ = fdxin^'n''R^,{xip)p^,p,, 20^'""' [xi) , (6.5) 



where 



dA 
— ( 
27r 



'IT 



5«^F,^^(0)Fr(An)} \p) 



(6.6) 



15 



In the case of H2 we get the additional factor (fci — xipY from the first-order term in Eq. 1 
This factor can be written as a derivative appHed to the exponential in (|2.4| ) if we introduce the 
tensor W\ 



a: 



(6.7) 



which subtracts the longitudinal component of the momentum-vector it is applied to. Now, us- 
ing partial integration in order to apply this derivative to the gluon operator and making use of 
Eq. ( |3l^ ) we get 



Ho 



dxldX2Rf^au{xiP, X2P) 



dX dry 
27r 27r 



giAxigi,?{x2-xi)^a 



IT 



{A'^MK'c{vn)A^a{^n)]\p), (6.8) 



where the partial derivative within the covariant derivative is due to the first-order term in the 
Taylor expansion (p. 7]), and the gluon operator within D^^ is due to the zeroth-order term in ( p.8[ ). 
For the derivation of (|6.8| ) we have also made use of the fact that, within the axial gauge, the 
a;-tensor applied to a gluon-operator yields the same operator: 



(6.C 



Projecting out the twist-four contribution with the help of Eqs. ( 5.11| ) and (5.13) and applying the 
method developed in section ^ in order to substitute the ^-operators by F-operators we end up 
with 



H2\r=A = ydxidx2 l^d'"'n''Rfj,aiy{xip,X2p) ]^d^iyip„i sG^'"'"' {xi,X2) + perm.j , (6. 



10) 



where 



3G'^'^^(xi,X2) 



UaTip f dX dry 



Xi X2 



__^,iAx,gi,(..-.0^^,(p|T {Ff{0)D:Uvn)Fr{Xn)] \p) . (6.11) 



For the permutation in ( |6.10| ) one has to keep Eqs. ( ^.11 ) and ( 5.13| ) in mind, i.e. indices of the 
projector which is applied to R/j^au and of the projector which is applied to ^G^ " ^ have to be 
permutated simultaneously. Similar to B.\^ (|6.5| ), 112 contains only contributions of the two-gluon 
operator. 

We expect that Eqs. (|6.5| ) and ( |6.10 ) are of more general validity, since the Ward identities we have 
used are valid to all orders and our ansatze for the color structure of iij^^ and R°'^xv> ( 3.15 ) and 
( 3.19| ), are of the most general nature. 



6.2 Four Gluons 
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6.2.1 Basics 



The situation is more complicated in the case of H3. Putting together ah terms which enter into 
and making use of Eqs. (|339| ), (|330D and (|33ll ) we get 



(6.12) 



where 



+ ^9f"''' RfMaxA^lP^ Xsp) sG'^'^^'^^ixi, X3)6{xi - X2) + 

+ ^9r'"^Rj,axAxiP,X2P,Xsp)3G'^^^'^{xi,X2)d{x2 - X3) + 

+ \ ^'"^ [^laXu (XIP, X2P, X3p) + 

+ Rla\uixiP,X2P,Xsp)] 2G^^^''(xi)5(xi - X2)6{x2 - X^)^ 



(6.13) 



^G^abdl{xi,X2) 



d^d^^iAx,^i,(..-.,)(p|T|^;^(0)ia^'K(r?n)^(^(An)]} 



( 

2tt 



iXx 



^p\T{A'^Mid^'id^'A^,{Xn)}\p). 



The derivatives in Eqs. ( 6.14 )-( 6.16 ) have to be understood as d'^ A{Xn) = dA{z)/dZf^r 
Due to our definition (3.25), we can spht up Eq. ( |6.12| ) into three parts 



(6.14) 
(6.15) 
(6.16) 

=An) etc. 

(6.17) 



where corresponds to those terms in ( 6.12| ) which contain the function R'^^xy Since we intend 
to write the expression for H3 in a form which is exphcitly gauge invariant (similar to ( |6.5D and 
(|6.10| )) we have to modify it in such a way that partial derivatives occur only within covariant 
derivatives. This can most easily be done with the help of the following relation: 



r''A^,{Tn)d'' = r'^A^,{Tn)B:,^ir^n)+gr'"^'r'^A^,{Tn)AUvn). 



ca'cd a\ 



aba' ra'cd a A / 



(6.18) 



Inserting ( |6.18 ) into ( |6.12| ) we realize that the results for the functions consist of terms with 
two-, three- and four-gluon operators. In the spin-averaged case the three-gluon operator does not 
contribute. The decomposition ( |6.17 ) is not unique, since the Ward identities could be expressed 
with the help of different combinations of the functions R'^^xu ™ -^l^- ( 3.29| ), ( |3.30 ) and ( 3.31D . We 
could, therefore, also write Eq. ( |6.12| ) with different combinations of these functions. By exploring 
symmetry relations for R'^fl'^\y{ki, ^2,^:3) we will show how to derive expressions for H3 which involve 
only one of the functions R and contain no three-gluon operator. 
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RlaxAki,k2,k3) = 


RlaxAkuk2,k3) - 


1 




Rlt^xv^h - ki,k2,k3) 


Rl^X.ik2-kuk2,k3) 


2 


Rlauxikl,k2,k2~k3) 


Rla^xif^l^k2,k2 - fcg) 


Rla^xikuk2,k2~k3) 


3 


RlxaAkuki^k2+k3,k3) 


Rlxauikukl-k2+k3,k3) 


Rixaui^l^^l -k2 + k3,k3) 


4 


Rlx,,uif^2 -ki,k3~ ki,k3) 


Rlx^Af^^ ~ki,k3- ki,k3) 


Rlx^Af^2 - ki,k3 - ki,k3) 


5 


Rlf^aui^s - k2,ki - k2 + k3, fca) 


Rxp^aui^s - k2,ki ~ k2 + k3, ks) 


Rxp^ai^i^S ~ ^2,^1 - fc2 + fcs, ks) 


6 


Riat.,.ik3 -k2,k3- ki,k3) 


RxatJ.u{k3 -k2,k3- ki,k3) 


Rxap.u{k3 -k2,k3- ki,k3) 


7 


RUxaif^l^>^l - k3,ki - fcz) 


RluXaif^l^>«l - k3,ki - fcz) 


Rl^Xaiki,ki " k3,ki - fca) 



TABLE 1: Making use of symmetry relations for ^J^^'Jif^Cfci, ^2, ^3) one finds relations between the three 
functions K'^^^^{ki,k2, ks). 



6.2.2 Symmetry Relations of the Four-Gluon Quark-Loop Expression 



The expression for the four-gluon quark-loop, R'^^\y{ki,k2,k^)-, has to be completely symmetrical 
under the simultaneous exchange of gluon- momenta, Lorentz- and color indices. This can be seen in 



two ways. First, within the path- integral formulation we have defined R^^^x^iki, k2,k^) in Eq. ( 3.12| ). 



Since the gluon is a boson the derivatives in ( |3.12| ) commute and the proclaimed symmetry follows 
immediately. Secondly, it can be seen diagrammatically. Clearly, R'^fl'^\y{ki,k2,k^) corresponds to 
the expression which is obtained if one takes into account all diagrams up to a specified order, i.e. 
in leading order one has to consider all permutations for the coupling of four t-channel gluons to 
quarks and anti-quarks. Consider, for example, the diagram in Fig. If we set ki = k2 — k[, 
^2 = ^2' ^3 = ^3 ^iid exchange /i ^ cr and a 6 we get an analytic expression which is the same 
as the one we get from the diagram in Fig. Since such a relation can be found for each diagram 
that enters into R^^^^y we have in general 



RtiuikiMM) 



R. 



bacd 
a^Xu 



{k2 - ki,k2,ki) 



(6.19) 



Next, with the help of the color structure (|3^ ) of Rf^^^ we can find relations between the three 
functions R'^^^^ if we realize that for a ^ h: d"'""^ d"'"^", d"^'^" d"'""^ but d"'''"^ d""'"^, as 
follows from the definition of d"'^'^'^, (3.24). With these relations we get 



RlaXuiki,k2,k3) = Rl^xuik2 



R 



Rl„Xy{ki,k2,k3 



,j,aXuih,k2,k3) 
3 

furXiy 



afiX. 
RlfiXu{k2 



ki, k2, k^) 
ki,k2, ks) 
ki,k2, fcs) 



3.20) 
3.21) 
3.22) 



Exchanging other gluon lines one gets a large number of relations. Those which are useful in the 
subsequent derivations and a few additional ones are given in Table 



6.2.3 Gauge Invariant Expressions 

With the help of the relations in Table |l] we can, now, write Eq. ( |6.12| ) as an explicitly gauge 
invariant expression. For example, using the relation in row one, column two of Table |l] we can, in 
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(a) (b) 

FIG. 4: Two sample diagrams which enter into They are identical if one sets fci = k'2 — k[, fc2 = 

— fcg and exchanges fi ^ a and a ^ b. 



the first line of ( 6.12| ), replace the function R^^xu ^a/iXu ^^'^ after an appropriate transformation 
of integration variables (xi X2 — xi, A r?) and an exchange of Lorentz- and color indices Rj^^xu 
can be factored out. Similarly the terms involving R^^^xv ™ (|6.12| ) can be brought into the desired 
form with the help of the expression in row three, column three of Table ||. We end up with 

i/3 = |dxidx2dx3 j ^^^e'^^ie'''(^^--'i)e-(-'3--2)4,,,(xip,X2P,X3p)a;^,o;^, x 

X (p|T{^^(0)[ff2|f^'^Af (rn)<M - iD^:,(rn)D^;,(r?n)]^(^(An)}b) , (6.23) 

where 

jabcd jadc' jc'bc _j_ ^ ^jafe^cd _j_ ^ac^bd _j_ ^ad^bc^ 

^ race' rc'bd ^mcc'ic'bd \ }_ f ^ab^cd ^acdid , eac^r6c^ i ^ / X'^^^'^'^ _j_ A'^'^A^'^ 4- A'^'^A^'^ 

^ race' rc'bd , n^abcd 

~ 2 



From Eq. (§J^) we realize that there is a term proportional RI^J^xu ™ ^'^^ ^^^^ -'-'^^ ( |3.25D 



this yields terms with color factors d"'^'^'^ and corresponding permutations of color indices. For the 



derivation of (3.23) one has to add and subtract terms which involve four yl-operators in order 
to get covariant derivatives. Combining the subtracted expressions with those in the first line of 
Eq. (6.12) the resulting color factor is no longer proportional to d"''"^'^. Instead one gets the term in 



Eq. (|6.2^ ). 



Again, with the help of ( |5.15| ) and (5.16) we can project out the twist-four contributions and with 



the help of the results in section |^ we can substitute the ^-operators by F-operators: 

^^31^=4 = ydxidx2dx3{c?^'^d'^^4, Xu i^iPi xsp) d^v'tr'A' ^ " {xi, X2, X3) + perm. I , 

(6.25) 
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where 



4G'l'^^''{xi,X2,X3) 



dX dr] dr 
27r 27r27r 



giAxigi?7(x2-a;i)gir(a:3-X2) ^ 



(6.26) 



-1 



XIX3{X2 - Xl){x3 - X2) 
1 1 



n^npn,ns{p\T{g''lf'^''Ff{0)F^^'yirn)F;;^ir^n)Fr{Xn)}\p) 



-n„n^(p|T{Ff (0)a;^,D^:(rn)u;^,D-;(7?n)Fr(An)}|i.) 



2 X1X3 

and the permutations in ( 6.25 ) have to be understood in the same sense as in ( 6.10| ) (cf. (|5.15| ), 
(|5.16| )). For H3 we get contributions of the two-gluon operator as well as of the four-gluon operator. 

The representation of H3 in a gauge invariant form as in Eq. ( |6.25| ) is not unique. For example, 
we can derive expressions where the quark- loop is represented by or instead of R^. These 
expressions can be obtained in two ways. One can either start with Eq. ( 6.12| ), make use of the 
relations in Table |l] and apply the same methods which led to ( |6.25| ). Or one starts with ( 6.25 ) 
and rewrites R^ in terms of R^ or R^ with the help of the expressions in Table |^. The results 
can be summarized as follows. If we exchange R^^^y^^{xip,X2P,xzp) by i?^^^^^(xip, X2P, xap) or 

R^^^^j^{xip, X2P, xsp) in ( |6.25 ) we also have to exchange aGi'^^'^ {xi, X2, x^) by aG^^'^^'^ {xi, X2, X3) or 

aG^'^^'^ {xi,X2,X3), resp. The latter functions can be written as 



AGi^''^''{xi,X2,Xs) 



dX dr] dr 
27r 27r27r' 



giAxigi?7(x2-a;i)giT(a:3-X2) ^ 



(6.27) 



X1X3{X2 - Xi)(x3 - X2) 



n, 



.n^n,ns{p\T{g'lt'^''F;^\0)F,^''{Tn)F^^{r]n)Fr{Xn)}\p) 



1 



n,r^^(p|T{F/(0)^,^D^:(rn).;;:,D^,(An)F^(^n)}b) 



2 {X2 - Xi)x3 



and 



aG'^''^''{xi,X2,xs) = j 



dX dr] dr 



,iAxi i?7(x2-xi) ir(x3-X2) 



-1 



X1X3(X2 - Xi)(x3 - X2) 
1 1 



27r 27r 27r 

ri^ripri,ns{p\T{gHf'''F-,\f))F^^{Tr^)Ff{r]ri)Fr{\n)}\p) 



(6.28) 



2 X1X3 



n,n^(p!T{F;'3(0)u;^,D^:(,?n)a;^,D^;(rn)Fr(An)}|p) 



where 

labcd 
In 



and 



labcd 
'3 



(6.29) 



adc' fc'bc 



^^adc'^c'bc ^ W ^ab^cd ^ gacgbd _ gadgbA ^ 1 / ^ab^cd ^ gacgbd ^ gadgbc 

2 2 V / 4 V 



_ jade' jc'bc _^ ^ ^ 5""^ 5^'^ -\- + J^'^j'"^ 



(6.30) 



1 O 1 

_ ^abc' rc'cd _jabc' Tc'cd , _ 

2 2 2 



_^ab^cd _|_ ^ac^bd _|_ ^ad^bc 



4 
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The color factors If'^'^'^ are related by permutations of color indices: 

labcd ^ labdc ^ lacbd _ (g g^) 

Note the different Lorentz- and color structure and the different dependence on integration variables 
of the two-gluon operators in (|6.26 ), ( 6.27 ) and (|6.28| ). 



6.2.4 Color Structure 

As we have seen, contains contributions of the two- and four-gluon operators. In this subsection 
we will explore the color factors Z?*'^'^, which enter into the contribution of the four-gluon operator 
and are given by Eqs. (|]2|), ( |639D and (|63o| ) depending on which representation of we are 



considering. If we have a look at Table we realize that there are also relations which transform 
one of the functions i?* into the same function BJ' and only the dependence on the momenta and the 
Lorentz-indices changes. For example, if we start with the four-gluon operator part of ( 6.25| ), make 



use of the relation in line seven, row one of Table ||, exchange the integration variables such that 
X2, — > X1 — X2, X2 X1—X3, X — > X—rj, r] —r] and r — > t—t] and make use of Lorentz-invariance we 



get an expression which has exactly the same form for the four-gluon operator part as in Eq. ( 6.26 ) 
but the color factor has to be exchanged by If'^'^'^, i.e. indices b and d are exchanged. Now, the color 



factor If""^ in ( |6^ has been written in three ways. There is the compact representation in the 



first line of ( |6.24D , and in the second line it has been split up into symmetric and anti-symmetric 
parts with respect of exchange of color indices a and c or 6 and d. The anti-symmetric part is given 
by the expression proportional f"''^'^ jcbd^ Since the expression we got after use of the relation in 
Table |^ must be the same as in ( |6.25D , we conclude that the term proportional f"''^'^ jcM -^^ (6.24) 



within (|]2D must give a vanishing contribution. In other words, the color factor that survives has 
the form 



We can repeat this argument with those expressions in which H'^ is expressed with the help of B? 
or i?^, ( 6.271) , ( |6.28| ). Again, the surviving color tensors have the form 



if'^'^ = M''^'^ , (6.33) 
if'"^ = 3d'"''"^ . (6.34) 



6.3 Discussion of the Results 

In the previous sub-sections we have presented our main results on the gluonic twist-four contri- 
butions, viz. Eqs. (I]!), (|6T0| ) and d^H). While the color structure of the expressions for Hi and 
H2 is valid beyond one loop^ we have some restrictions in the case of H^. For a four-gluon system 
in the color singlet, there are five linearly independent states in the case of even C-parity. For 
example, one can expand the color structure with the help of the five states 1, 8a, 83, 10 -|- 10 and 

*The Ward identities derived in section ^ are valid only for IPI vertex functions. Therefore, contributions of 
diagrams which involve three- and four-gluon vertices like those in Fig. ^ have to be incorporated separately into the 
expressions for Hi through H3. 
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27. However, a closer inspection of our ansatz Eq. ( |3.25| ) for the color structure of the four-gluon 
quark-loop shows that the state 10 -|- 10 does not occur and there are, in addition to the state 8a, 
only two specific linear combinations of the three states 1, 8s, and 27. In an analysis valid to all 
orders one should expect that all color states contribute. 

Let us, now, compare our results with those obtained by EFP As mentioned in the introduction, 
the computations in Q can be divided into two steps. In the first step one considers the relevant 
diagrams and performs the factorization of the coefficient functions and matrix elements of the 
relevant operators. This leads to the operators in Eqs. ( |1.2|)-(|1.5|) . They can be compared with the 
two- and four-gluon operators in Eqs. (6^), ( |6.11] ) and ( |6.26D , which involve the same number of 
covariant derivatives as the corresponding quark operators. 

In the second step, EFP used equations of motion in order to express the contributions of the 
two-quark operators (1.2)- ([L^ ) in terms of linearly independent quantities. Within the transverse 
basis, the resulting twist-four contributions to Ft and Fl are given in Eqs. ( |1.8| ) and (^) plus 
terms which involve the four-quark operator (|1.5|). In our case this second step has not yet been 
performed. It cannot be done without an explicit computation of the quark-loops R. 



On the other hand, one should try to find the relation between our results and the operators ( 1.10 ) 
and ( |1.11| ) occuring in the OPE. For this, one has, again, to compute the expressions for the quark- 
loops R. As has been described in Q, the connection to the operators ( |1.1C ) and ( |l.ll| ) can be 
established for moments of contributions to the structure functions with definite signature only, i.e. 
the expressions for ^G^'^" and j^G'^'^'^^ in ( |6.1l| ) and ( |6.26 ) have to be separated into symmetric and 
anti-symmetric parts. Note that the Lorentz-indices of the operators in ( |6.11 ) and ( 5.26| ) are not 
symmetrized. The permutations involved in Eqs. ( |6.10D and ( |6.25 ) imply the simultaneous exchange 
of Lorentz-indices of the quark-loop expressions and of the operators (cf. ( 5.11| ) and ( ^.15D ). Taking 
moments of Eqs. (|6.5|), ( |6.10| ) and ( |6.25 ) will then result in an expansion of the operators in ( |6.6[ ), 
(|t|) and (|6l2^ ) about the origin, so that finally the operators (|1.10D and (|1.11[) will be recovered. 
In summary, the approach used here is rather intuitive since one starts with certain diagrams which 
yield properties of the coefficient functions which, in turn, imply restrictions on the contributing 
operators. Within the OPE, on the other hand, one starts with a certain class of operators and 
has to find the corresponding coefficient functions. 

As mentioned before, the computation of the quark-loop expressions R will yield terms which are 
proportional to logarithms in and terms which are just constants. The terms which are propor- 
tional to logarithms have, actually, to be interpreted as part of the evolution of the contributions in 
Eqs. (p^ 



and (1.7). They are due to configurations with collinear momenta of quarks and gluons 
and can be interpreted as contributions where one of the gluons splits up into a quark-antiquark 
pair. Non-logarithmic terms, on the other hand, correspond to the NLO twist-four gluon coefficient 
functions. For the computation of the quark-loop expressions R it will, therefore, be necessary to 
re-identify the diagrams considered by EFP. 



7 Conclusions 

For the description of deep inelastic structure functions in the low-Q^ regime, the determination 
of the higher- twist contributions remains one of the outstanding theoretical challenges. At least in 
the small-XB domain, some of these higher-twist contributions are known to increase faster with 
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decreasing xb than leading-twist terms. In addition, their l/Q^-suppression is weakened due to 
evolution which is stronger than in the leading-twist case |^. In this paper we have presented results 
on gluonic twist-four operators, which are expected to give the main twist-four contributions at 
small-XB- For the derivation of these results we have extended the method developed by EFP Q. 
Based on the diagrams in Fig. |^ and working in the axial gauge we could perform the factorization 
of the upper quark-loop parts and the lower blobs. The latter are expressed as matrix-elements 
of gluonic operators between proton states. This factorization procedure requires the use of Ward 
identities, which have been derived in section ^. The axial gauge proves to be very convenient for 
two reasons. First, gluonic j4-operators occurring between proton states can easily be expressed in 
terms of F-operators which is necessary in order to obtain gauge invariant expressions (section 
Secondly, diagrams with more than four t-channel gluons lead to contributions of, at least, twist-six. 
This statement is proved in section ||. Working in a covariant gauge, one would have to take into 
account contributions of diagrams with arbitrarily many t-channel gluons. 

After the factorization had been performed we ended up with three contributions. The two- and 
three-gluon quark-loops are convoluted with matrix elements of the two-gluon operator, Eqs. ( |6.5D 
and ( |6.1[1|) , while the four-gluon quark- loop is convoluted with both, the two- and the four-gluon 
operator, Eq. ( |6.25| ). In order to derive the latter result, we had to make intensive use of symmetry 
relations for the expression of the four-gluon quark-loop. These symmetry relations also led to 
restrictions on the color structure of the four-gluon operator contribution. As the next step we 
have to define a minimal set of operators in terms of which we can express the operators in ( |6.5D , 
( 6.10| ) and ( |6.25| ). This requires the introduction of signature of the two- and four-gluon t-channel 



states, i.e. a further step is symmetrization of the color and Lorentz indices. It is only after this 
decomposition into contributions with definite signature that we can return to ( 1.10| ) and (1.11) 
and determine contributions to the structure functions. 

It is worthwhile to compare the results presented here with some earlier work. The traditional 
approach to small-xe physics is the high-energy (small-xe) approximation, which, in leading order, 
leads to the BFKL equation The BFKL equation is usually described by a ladder diagram of 
two reggeized t-channel gluons. It sums all leading logarithms in xb and contains contributions of all 
twists. Contributions of diagrams with three and four reggeized t-channel gluons have been studied 
in the generalized leading logarithmic approximation in |19,2^|^. The results of these calculations 



have been used |p to compute the twist-four part in the double- logarithmic approximation (DLA), 
and it is interesting to compare them with the results of the present paper. 

First, in Q it was shown that the three-gluon amplitude can be written as a sum of expressions for 
the two-gluon (BFKL) amplitude. Since this phenomenon generalizes the reggeization of the gluon, 
it has been termed reggeization, as well. Color structure and reggeization are independent of the 
twist. The fact that the three-gluon amplitude can be written as a sum of two-gluon amplitudes is 
therefore also true for the twist-four part. This result can be confronted with our new result that 
the three-gluon quark-loop couples only to the two-gluon operator ( |6.10D . Next, the four-gluon 
amplitude studied in |20| splits up into two terms: there is, again, a reggeizing part, which can 



be expressed with the help of expressions for the two-gluon amplitude, and a second part which is 



irreducible with respect to reggeization. We are, now, tempted to draw a relation with Eq. (6.25) 
which also contains contributions from two-gluon and four-gluon operators. In particular, one 
might associate the second matrix element in ( |6.25 ) with the reggeizing part of |20|, whereas the 



Recently, C.Ewerz |2l| has performed an analysis of the five- and six-gluon amplitudes. 
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first term in ( |6.25|) looks analogous to the irreducible four-gluon amplitude. Here, however, a word 
of caution is in place. First, in [^] it is the reggeizing part, which contains color factors like d"''^'^'^^ 
while the irreducible part is symmetric under simultaneous exchange of color and momentum. This 
is in contrast to the color structure of two- and four-gluon operators in ( |6.25| ). Secondly, the 
calculations use different approximations. Ref. |^0| studies the leading-log (I/xb) approximation. 
This has lead to the result that the lowest-order diagram with four i-channel gluons belongs to the 
reggeizing part, which suggest that there is no direct coupling of the four-gluon operator to the 
quark loop. In the present paper we have taken a different route: our analyis of corrections 
has not been restricted to the small-XB approximation, and we find contributions of the two- 
and four-gluon operator already at leading-order. This can be considered as a hint that one will 
encounter simplifications if one uses explicit expressions for the four-gluon quark loop in those 
terms which are coupled to the four-gluon operator in Eq. ( |6.25D and then takes the limit xb — »■ 0. 
Agreement between the results of and the present approach holds only in DLA; to go beyond 
this approximation is one of the main motivations of the present paper. 

In summary, we have performed the first step towards a consistent treatment of gluonic twist- 
four contributions at small-XB- What remains to be done is, first of all, the computation of the 
expressions for the quark-loops. While this task might be straightforward, one has to keep in mind 
that the expressions will contain contributions proportional to \ogQ^ and, therefore, contribute to 
leading-order quark-coefficient functions convoluted by Altarelli-Parisi splitting functions as well 
as to next-to-leading order gluon coefficient functions. A major task will be the computation of 
anomalous dimensions of the gluon operators. Since several gluonic operators contribute to twist- 
four, one will have to handle also the problem of operator mixing. First results are contained 



in [g] and [15|, but a complete LO analysis does not exist. Moreover, although we have argued 
that at small-XB gluon-operators are expected to dominate, the role of quark-operators has not 
been clarified. For a fully consistent treatment of twist-four contributions one should also take 
contributions of quark-operators into account. In such an analysis one will encounter the problem 
of operator-mixing, as well. We feel that calculations along these lines are urgently needed: a more 
complete understanding of deep inelastic ep scattering at low and small xb cannot be reached 
unless we know the influence of higher twist. 

Acknowledgements: We wish to thank L.N. Lipatov and R.K. Ellis for stimulating and very 
helpful discussions. 
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